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ABSTRACT 

Young Sun-like stars at the beginning of the pre-main sequence (PMS) evolution are 
surrounded by accretion discs and remnant protostellar envelopes. Photometric and 
spectroscopic variations of these stars are driven by interactions of the star with the 
disc. Time scales and wavelength dependence of the variability carry information on 
the physical mechanisms behind these interactions. We conducted multi-epoch, multi- 
wavelength study of PV Cep, a strongly variable, accreting PMS star. By combining 
our own observations from 2004-2010 with archival and literature data, we show that 
PV Cep started a spectacular fading in 2005, reaching an /c-band amplitude of 4 mag. 
Analysis of variation of the optical and infrared fluxes, colour indices, and emission line 
fluxes suggests that the photometric decline in 2005-2009 resulted from an interplay 
between variable accretion and circumstellar extinction: since the central luminosity 
of the system is dominated by accretion, a modest drop in the accretion rate could 
induce the drastic restructuring of the inner disc. Dust condensation in the inner disc 
region might have resulted in the enhancement of the circumstellar extinction. 
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1 INTRODUCTION 

Eruptive young stars, defined by their episodically increased 
optical and infrared brightness, form a small subclass of 
Sun-like PMS stars. Traditionally they are divided into two 
groups: FU Orionis type stars (FUors) exhibit an initial 
brightening of 5 mag during several months or years, fol- 
lowed by a fading phase of up to several decades or a cen- 
tury. The other group, EX Lupi-type stars (EXors) are char- 
acterized by relatively short, recurrent outbursts which last 
some weeks to months, and the time between the eruptions 
ranges from months to years. In both classes the outbursts 
are believed to be powered by enhanced accretion from the 
circumstellar disc on to the star. In the standard picture the 
inward spiralling material piles up close to the inner edge of 
the accretion disk, and falls on to the stellar s urface as a 
resul t of a gravitational and thermal instability l|Zhu et al.l 
I2OO9I ). 

* E-mail: kun@konkoly.hu 



In 2003, the outburst of the low-mass young stel- 
lar object V1647 Ori triggered a large number of obser- 
vations from X-ray to far-infrared wavelengths. Based on 
these ground-based and space-born measurements it was 
suggested that this object fits neither into the FUor nor 
EXor class. It was outlined that the outburst - unlike in 
other young eruptive stars where enhanced accretion is the 
main physical driver behind the brightening - consisted of 
the combination of two effects of comparable amplitude. 
They are an intrinsic brightening related to the appear- 
ance of a new, accretion-related hot component in the sys- 
tem, and a dust-clearing event which reduced the extinc- 
tion along the line-of-sight (e .g. lAcosta-Pulido et al. I l2007l : 
lAspin. Beck fc Reipurthll2008l ). The simultaneity of the ac- 
cretion and extinction changes suggests that they might 
be physically linked, and the changing accretion luminosity 
causes changes in the inner disc structure. Anot her example 
for su ch combined processes may be V1184 Tau IXi ves et al.l 
I1997I ) whose optical and infrared photometric observations 
have shown that the large photometric decline was associ- 



ated with variation in the inner disc structure l|Grinin et alj 
|2009|) . That the variable extinction may have an important 
role in the brightening refines the picture of the FUor/EXor 
phenomenon, moreover such rearrangements of the inner 
disc structure have a potentially high importance for the 
evolution of the terrestrial zone of circumstellar disk. 

The target of our multi-wavelength variability study, 
PV Cep is an eruptive young star, wh ich was origin ally 
classified as member of the EXor class bv lHerbid (Il989h. on 
the b asis of its large outburst in 1977-1979 ( Cohen et alJ 
Currently the star undergoes significant brightness 
changes which we have been monitoring for years. Our 
first observations indicated that PV Cep might be another 
example for accretion-induced structural changes, and this 
motivated our further detailed study of the object. The 
role of variable circumstellar extinction was also apparent 
in the dramatic brightness decline obser ved in the near 
infrar ed between 2008 April and June by iLorenzetti et al.l 
l|2009l ). Several further observed properties distinguish PV 
Cep from EXors, and indicate that the object may be 
more similar in nature to V1647 Ori. Such properties are 
e.g. the time-scale of the outburst (~ 2 years), the A-type 
absorptio n spectrum observ ed at the beginning of the 
outburst (jCohen et al.lll98l | ). the Class I spectral energy 
distribution (SED) ( Connellev et al.l |2008|), the associate d 
high-mass circumstellar disc ( 0.76 Mf7 ^, Hamidouchc '2010*), 
varia b le refl e ction nebula (|Gvulbudagliian fc Magaki an, 
19771: ICohenI \l98d). m o lecular outflow iKvreaultl Il984 
Arce k GoodmanI booA iHamidouch^ |201G|). and parsec- 



scale optical iet "(|G6mez. Whitnev &: KenvorJ 1 19971 : 
iReipurth. Ballv fc Devin^ I1997': 'Nec kel et al.lll987l v 

We carried out a comprehensive optical and infrared 
photometric and optical spectroscopic monitoring of P V Cep 
during the low-brightness state in 2008-2010, and we have 
a few additional measurements from 2004, 2005, and 2006. 
In this paper we combine these data with archival and new 
near- and mid-infrared data, obtained in 2004-2010, to bet- 
ter understand the nature of PV Cep and the nature of 
physical processes which affect the inner disk structure and 
evolution. 



2 OBSERVATIONS 
2.1 Photometry 

Photometric observations in the VRcIc bands spanning 
the time interval between 2004 September 22 and 2010 
November 23 were performed with six instruments on 
flve telescopes. Most of the data were obtained with the 
60/90/180 cm Schmidt telescope of the Konkoly Observa- 
tory, equipped with a Photometries AT 200 camera and 
with the 1-m RCC telescope of the Konkoly Observatory, 
equipped with a Princeton Instruments VersArray:1300B 
camera. PV Cep was also observed with the IAC-80 tele- 
scope of the Teide Observatory (Spain) between 2009 Octo- 
ber 25 and November 7. More detailed description of these 
instruments is given in lAcosta-PuUdo et al.l (|2007l ). In 2008 
August and 2009 October we obtained VRI images with the 
CAFOS instrument installed on the 2.2 m telescope of the 
Calar Alto Observatory (Spain). Exposure times were be- 
tween 90 and 180 s per image. The Ic = 11.70 ±0.04 magni- 
tude on 2004 September 22 was derived from a CCD-image 
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Figure 1. V, Rc, Ic, J, H, Ks, and IRAC [3.6] and [4.8] light 
curves of PV Cep between 2004 May and 2010 September. V 
magnitudes in 2008 February were taken from the INTEGRAL 
OMC Archive. N ear-in frare d data before 20 9 are taken from 
ILorenzetti et al. I 1I2OO9I) and IConnellev et al. I 1I2OO8I '). The light 
grey dashed lines, spanning the long time-gaps without data, were 
drawn to clearly identify the coherent data sets. Arrows indi- 
cate the dates of our spectroscopic observations. Roman numbers 
and shading indicate the five characteristic segments of the Ught 



obtained with the 2.6-m telescope of B yurakan Observa- 
tory, Armenia, and shown in the paper of ' Movsessian et al.l 
(2008). The image of PV Cep was saturated, and its magni- 
tude was derived by fitting the PSF to the unsaturated wing 
of the stellar image. On 2010 April 8 we used the Electro 
Muhiplying (EM) Andor Technology iXon"^^ -I- 888 cam- 
era at the RCC telescope of the Konkoly Observatory. We 
obtained a series of 100 exposures in the Ic band with expo- 
sure times of 20 s per image in EM mode (the nominal EM 
gain was 50). The image processing was made with the soft- 
ware tools written by Pal (2009). All the other images were 
reduced in IRAF. In order to transform instrumental mag- 
nitudes into the standard system, we calibrated 16 stars in 
the field of view of the 1-m telescope (7' x 7') in VRcIc 
bands. Calibration was made during six photometric nights. 
Standard stars in NGC 7790, published bv iStetson (i|2000l ), 
were used as reference. Equatorial coordinates and derived 
VRcIc magnitudes with uncertainties of the comparison 
stars are listed in Table [T] The results of the photometry 
for PV Cep are presented in Table [2] and plotted in Fig. [1] 
The photometric errors were derived from quadratic sums of 
the formal errors of the instrumental magnitudes and those 
of the coefficients of the transformation equations. The tele- 
scope used for each observation is shown in the last column 
of Tabled 



2.2 Near-infrared observations 

Near-infrared, JHKs measurements of PV Cep were ob- 
tained using the NIC-FPS camera on the ARC 3.5-m Tele- 
scope at Apache Point Observatory on 2009 October 12 



Table 1. Comparison stars* 



N 


RA{2000) 


D(2000) 


V{AV) 


Rc (ARc) 


Ic (A/c) 




h m s 


d ' " 


[mag] 


[mag] 


[mag] 


CI 


20 46 28.89 


+67 59 05.4 


14.072 (0.012) 


13.595(0.012) 


13.200(0.008) 


C2 


20 46 26.50 


+67 58 11.8 


14.608 (0.005) 


13.948 (0.004) 


13.495(0.005) 


C3 


20 46 10.44 


+67 55 44.2 


15.134 (0.005) 


14.205 (0.004) 


13.637(0.004) 



*The full Table is available as Supporting Information with the online version of this article. 
Table 2. Results of the photometry of PV Cep** 



JD 


Date 


ViAV) 


Rc (ARc) 


Ic (A/c) 


Telescope 


(2454000+) 


yyyymmdd 


[mag] 


[mag] 


[mag] 




-329 


20051027 


14.58(0.16) 


13.28 (0.06) 


11.78 (0.06) 


RCC 


-327 


20051029 


15.08(0.04) 


13.39 (0.04) 


11.90 (0.11) 


RCC 


-326 


20051030 


15.43(0.06) 


13.59 (0.04) 


11.98 (0.06) 


RCC 


-325 


20051031 


15.14(0.07) 


13.36 (0.01) 


11.84 (0.16) 


RCC 



**The full Table is available as Supporting Information with the online version of this article. 



(JD 2455116.5). Dark and domeflat images were also ob- 
tained, and four-point dithering was applied to ensure 
proper sky subtraction. To avoid saturation, the Ks images 
were slightly defocussed, and, to exclude a nearby bright 
star from the field of view during the H and Ks expo- 
sures, the camera was rotated. The images were reduced 
in IRAF. Magnitudes were derived by aperture photometry, 
and were transformed into the standard s ystem by compar - 
ing them with the 2MASS magnitudes jcutri et al.l |2003| ) 
of four stars in the image field. These stars are as fol- 
lows: 2MASS 20462889+6759054 (CI in Table [T]), 2MASS 
20462650+6758118 (C2), 2MASS 20461044+6755442 (C3), 
and 2MASS 20461128+6800012. The derived magnitudes 
are as follows: J= 11.66 ± 0.04, H= 9.42 ± 0.04, Ks = 
8.00 ±0.10. 



2.3 Spectroscopy 

Intermediate resolution spectra of PV Cep were obtained on 
2008 August 26 and 29, 2009 October 9, 11, and 14, using 
the CAFOS instrument on the 2.2-m telescope of the Calar 
Alto Observatory. The R-lOO grism covered the 5800-9000 A 
wavelength range. The spectral resolution of CAFOS, using 
a 1.5-arcsec slit, was A/AA « 3500 at A = 6600 A. The spec- 
trum of a He-Ne-Rb lamp was regularly observed for wave- 
length calibration. Broadband VRcIc photometric images 
were taken immediately before the spectroscopic exposures 
for flux calibration. We reduced and analysed the spectra 
using standard IRAF routines. A further R-lOO spectrum, 
obtain ed on 2004 August 7, and shown in fig. 2 o fH un et al.l 
l|2009l ), is also available. 



2.4 Spitzer data 

We observed PV Cep using the Spitzer Space Telescope be- 
tween 2004 October and 2005 August with the IRAC, IRS 
and MIPS instruments. In addition, we used further archival 



data from IRAC and MIPS. Moreover, we conducted a mon- 
itoring with the IRAC instrument during the post-Helium 
phase in 2009 September and 2010 January. A log of all 
Spitzer observations of PV Cep is presented in Table |31 
The data reduction procedures are described in Appendix A, 
available as Supporting Information with the online version 
of this article. The results of Spitzer photometry are shown 
in Table H 



3 RESULTS 

3.1 Photometric variations 

Light curves of PV Cep in 2004-2010. The light curves 
of PV Cep in the V, Rc, Ic, J, H, K^, and IRAC 3.6 
and 4.5 /im bands between 2004 May and 2010 November 
are shown in Fig. [1] Based on the shape of the Ic curve 
we divided the covered period into five segments, indicated 
by shading and Roman numbers. The first segment, there- 
inafter referred to as the bright state, is the interval between 

2004 July and 2005 October. Segment II, the fading phase, 
is the period between 2005 October and 2007 October, when 
the star faded significantly in each photometric band. Seg- 
ment III, the transtent peak, covers the brightening between 
2007 October and 2008 February, and the subsequent sharp 
decline until 2008 June. Segment IV is the low-brightness 
state between 2008 June and 2009 August, when the op- 
tical brightness slowly decreased further, reaching a mini- 
mum of Ic = 15.86 on 2009 August 18. After this date, in 
segment V, the optical and infrared fiuxes, except the Ks- 
band fiux, started rising, and large, short-term brightness 
fiuctuations appeared in the optical light curves. We refer 
to this interval as the rising phase. The total amplitude of 
the Ic light curve is about 4 magnitudes over the period 

2005 November and 2009 September. 

Colour— magnitude diagrams. In order to get insight 
into the origin of the observed variations we plotted colour- 



Table 3. Log of Spitzer observations of PV Cep 



Instrument 


Wavelength 


Date 


AOR 


Program 


IRAC 


3.6, 4.5, 5.7, S.OAtm 


2004 Oct 29 


11571712 


3716 


IRS 


5.2 - 37.2 /^m 


2004 Oct 23 


12287488 


3716 


MIPS SED 


55 - 95 /im 


2005 Aug 03 


11571968 


3716 


IRAC 


3.6, 4.5, 5.7, S.OAtm 


2006 Nov 26 


18955008 


30760 


MIPS 


24, rOfim 


2007 Feb 26 


19962624, 19962880 


30574 


IRAC post-He 


3.6, 4.5 /im 


2009 Sep - 2010 Jan 


35591168 - 35678720 


60167 



Table 4. Spitzer photometry for PV Cep. AU fluxes are color-corrected and given in Jyt. 



Date 




F4.5(jm 




Fg.OMm 




20041029 


4.47(0.09) 


6.98 (0.14) 


8.34(0.17) 


11.85 (0.23) 




20061126 


1.84(0.04) 


3.99 (0.08) 


5.45(0.11) 


8.16(0.17) 




20070226 










27.12(1.04) 35.22(2.13) 


20090916 


3.54(0.03) 


5.97 (0.02) 









tThe full Table is available as Supporting Information with the online version of this article. 



magnitude diagrams. Figure [5] (left) shows that the fading 
in 2005-2006 was nearly grey in the Rc and Ic bands. Its 
reason may be that in minimum, whatever is the origin of 
the dechne, the star itself is too faint to be detected, only 
the hght scattered from the disc atmosphere, thus bluer in 
colour, can be observed. The transient peak shows different 
colour variation: the star became redder and brighter and 
then bluer and fainter. This colour behaviour is character- 
istic of UX Orionis (UXor) type varia bles close to the hght 
curve minima (e.g. iBibo fc Th^ll99ll ). and is attributed to 
the increasing proportion of scattered hght when the star is 
obscured by a circumstellar dust clump. This diagram sug- 
gests that the optical source in the dim phases (segments II 
and IV of the hght curve) was the starlight scattered from 
the disc atmosphere. During the low-brightness phase (filled 
circles) PV Cep was redder when fainter: this behaviour can 
result from either enhanced extinction or vanishing hot spots 
on the stellar surface due to the decreasing accretion rate. In 
the rising phase the star turned brighter and redder, indicat- 
ing that instead of the scattered hght direct hght from the 
central star could be detected. The data obtained in 2010 
indicate that the star completed a full loop in this diagram 
from 2008 April to 2010 April. 

The near-infrared data in Fig. [2] (right) show that the 
fading in 2005-2006 was accompanied by a slight decrease 
of the H ~ Ka colour index, indicative of decreasing emis- 
sion from the dust sublimation zone of the disc, whereas an 
extinction change of Ay ~ 7 mag can account for the pho- 
tometric variations around the transient peak in 2008. The 
different colour behaviour suggests the different nature of 
the two brightness drops. The strong decrease of the H — Ka 
colour index at the bottom of the diagram, measured in 2009 
October, clearly suggests the fall of the inner disc emission 
during the low-brightness phase. 

Variations in the mid- and far infrared The Spitzer 
data obtained in the bright, fading, and rising phases of the 
light curve (Table [4]) allow us to inspect the wavelength de- 
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Figure 2. Left: Rq magnitude as a function of the i?c — Ic 
colour index. Consecutive data points are connected by dotted 
lines. The arrow indicates the colour dependence of the inter- 
stellar extinction. The typical error bars are shown in the upper 
left corner. Right: Ka magnitude as a function of the H — Ka 
colour index based on all data obtained betw een 2004 and 2010 
llConnellev et al.ll2008l: iLorenzetti et al.ll2009l) . Consecutive data 
points are connected. 



pendence of variations in the mid- and far-infrared regions. 
We plotted in Fig. |3] the flux ratios between the different 
segments of the hght curve as a function of wavelength. In 
addition to our own measureme nts, th e bright state data also 
include H and Ka data from i Connelley ct al. (2008), and the 
fading p hase data set contains JHKa data measured on 2007 
May 13 (ILorenzetti et aLll2009i ). as well as the AKARI IRC 
ijishihara et al.1 \2Q1^ fluxes at 9 and 18 um, and AKARI 
FIS catalog data (|Yamamura et all [2OI0I ) at 65, 90, 140, 
and 160 /im, each obtained between 2006 May 8 and 2007 
August 26. The wavelength dependence of variations is ob- 
viously different from that of the dust extinction: the excess 
emission in the bright state with respect to fading phase 
had an optical, a near- and mid-infrared, and a far-infrared 
component, indicative of changing emission from both the 
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Figure 3. Flux changes of PV Cep between the different segments 
of the light curve, as a function of the wavelength. The dotted 
line shows the flux ratios corresponding to an increase of the 
extinction by Ay = 5 mag llCardelli et al.lll989l. asuming Ry = 
3.1). 
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central star and various parts of the disc. The optical and 
near-infrared excess components are separated by a mini- 
mum in the H band. The wavelength interval 10-65 fim was 
unaffected by the fading, while the system also dimmed in 
the 70-90 /im region. The Fa (2006) /Fa (2009) curve shows 
that both the star and the inner disc continued fading after 
the transient peak. 

3.2 Variations in the nebula RNO 125 

Variations in the nebula RNO 125. RNO 125 was fan- 
shap ed and bright in the bright state (|Movsessian et all 
|2008| ). and the bright head of the cometary nebula was cen- 
tered on the star (Fig. 2] upper left). A dark band appeared 
between the star and the nebula during the fading phase 
(Fig. 31 upper right). The nebula quickly faded and exhib- 
ited a strongly changing shape after the transient peak. Sim- 
ilar phenomenon was observed in 1977-1979 (see fig. 1 of 
ICohen et al] |l98l"). The nebula was barely visible in our im- 
ages during the low-brightness phase and at the beginning of 
the rising phase (Fig. |4l lower left). The head of the nebula 
started brightening in 2010 (Fig. 2] lower right). 

The variable shape of the nebula suggests that the illu- 
mination changed not only due to the fading of the star, but 
also due to the changing geometry of the dust distribution 
close to the star, resulting in rapidly changing shadows on 
the matter inside the outflow cavity. The rebrightening of 
the nebula started close to the star: those parts appeared 
first which dimmed when the star started declining. 



Figure 4. I-band images of the nebula RNO 125 at different 
dates. The size of each image is 1' X 1'. North is up, and east is 
to the right. 

widths and fluxes of the strongest emission lines in the op- 
tical spectra of PV Cep, obtained in the bright, low-state, 
and rising phases of the hght curve. The fluxes were de- 
termined using the Rc and Ic magnitudes, measured si- 
multaneously with the spectroscopic observations, a nd the 
fluxes corresponding to zero magnitude tabulated bv lGlasd 
(|l999h . Since only Ic magnitude is available for 2004, we es- 
timated the Rc magnitude in 2004 using the average colour 
index Tic — Ic = 1.50, measured at the end of the bright 
state in 2005. The variations of the emission lines can be 
summarized as follows. (1) The equivalent width of the Ha 
line remained nearly unchanged during the photometric de- 
cline, that is, the flux of the Ha emission and the contin- 
uum flux underwent equal fading, suggesting that the star 
and the Ha emitting region were attenuated by a circum- 
stellar dust structure. (2) The fluxes of the forbidden lines 
decreased only slightly during the decline of the stellar flux 
(the equivalent widths strongly increased), indicating that 
the formation regions of these lines were unaffected by the 
changes leading to the fading. (3) The decline of the Call 
line fluxes was signiflcantly larger than that of the under- 
lying continuum: their equivalent widths were halved be- 
tween 2004 and 2008. Since these line s are accretion tracers 
(|Muzerolle. Hartmann fc Calvetlll998l ). the accretion rate of 
PV Cep dechned between 2004 August and 2008 August. 



3.3 Variations in the optical spectrum 

No photospheric absorption lines can be identifled in the 
spectra. With the exception of Ha, the strong emission lines 
are symmetric at our resolution. The forbidden lines are sys- 
tematically blueshifted with respect to the permitted lines 
by some 100-200 km s~^, indicating their origin from the 
approaching lobe of a wind. Table [S] shows the equivalent 



3.4 Physical properties of PV Cep 

Mass and luminosity. Given that the amplitude of the 
photometric variations was smallest in the //band (Fig. [S]), 
we estimated the luminosity of PV Cep from the H mag- 
nitude measured in the low-brightness state on 2008 June 
18, assuming that the total flux in the H band originated 
from the pho tosphere. We adopte d a distance of 325 pc 
for PV Cep (|Straizvs et al.1 Il992l ). an effective tempera- 



Table 5. Equivalent widths and fluxes of emission lines 



Id. EW (A) 

2004 2008 2009 2010 



F(10-" Wm-2) 
2004 2008 2009 2010 



Ha 


104 


132 


100 


90 


99 


6.6 


3.8 


9.3 


Call 8498 


77 


42 


62 


39.5 


168 


4.1 


6.3 


11.0 


Call 8542 


74 


39 


59 


37.4 


161 


3.8 


5.9 


10.4 


Call 8662 


64 


32 


51 


30.0 


139 


3.1 


5.1 


8.3 


[OI] 6300 


8.1 


58 


82 


37.0 


4.8 


2.9 


3.1 


3.8 


[OI] 6363 


1.7 


20 


27 


13 


1.6 


1.0 


1.0 


1.3 


[SII] 6717 


1.2 


4.6 


9.2 


3.3 


1.1 


0.23 


0.35 


0.34 


[Sll] 6731 


1.6 


8.6 


15.6 


5.9 


1.5 


0.43 


0.59 


0.61 



ture Ted — 5500 K, corresponding to th e spect ral type 
G8-K0 found by iMagakian fc MovsessianI (|200ll fl. a to- 
tal extinction Av = 12.0 mag, obtained by deredden ing 
the near-infrared colour indices (|Lorenzetti et al ]|2009^ on 
to the T Tauri locus of the J — H vs. H — Ks diagram 
l|Mever. Calvet fc Hillenbrandlll997l'). and applied the bolo- 



metric correction tabulated by Hartigan. Strom fc StromI 
and unreddened colour index V — H = 1.69 for the 
spectral type G8 (|Kenvon fc Hartmann|[l995l ). We obtained 
Lt, « 17 Lq, and a stellar radius J?« « 4.0 Rq. Plotting the 
re sults on the HRD, toge ther with the evolutionary models 
of iPalla fc Stahlej l| 19991 ), results in a mass of 2.4 Mq and 
age of < 1 Myr. 

Disc accretion rates. We estimated the low-brightness- 
state disc accretion rate of PV Cep from the luminosity of 
the Call 8542 line, measured in 2008 August (Table O, 
using the empirical formula of .Dahmi ((2008i 'l . We obtain 
« 41 Lq, and Mi°^ « 2.6 x 10"^ Mg/yr. The corre- 
sponding values for the bright state, estimated from the flux 
change of the Ca ll 8542 line, are about twice higher. The re- 
sults show that the accretion rate of PV Cep is significantly 
hig her than typical PMS accretion rat es (~ 10~*MQ/yr, 
e.g. iGuUbring etHI Il998l : iDahml |2008| ) both in the high- 
and low-brightness states, and the accretion luminosity is 
higher than the photospheric luminosity. 



4 POSSIBLE NATURE OF THE VARIATIONS 

The decreased fluxes of the Call emtssion lines of PV Cep 
between 2004 and 2008 suggest that the fading was associ- 
ated with decreasing accretion rate. The factor of two drop 
of the accretion luminosity, derived from the Call 8542 lines, 
however, can account for only < 1 mag variation in the op- 
tical region. To explain the whole amplitude of ~ 4 mag, 
enhanced circumstellar extinction also has to be invoked. 
The excess extinction could arise from the changes in the in- 
ner disc structure. Since Lace > L,, the modest drop of the 



^ Th e widely accepted spectral type of PV Cep, A5 l lThe et al.l 
[l99i), is based on the Balmer absorption lines, ob served in the 
spect rum at the beginning of the outburst in 1977 l lCohen et al.l 
Il98 j) . No spectrum, obtained after the outburst, exhibited 
Balmer absorption lines, which thus might have been originated 
from the hot inner disc in early outb urst phases, even as in the 
case of V1057 Cyg and V1647 Ori (cf. IWelinlll97ll: iBriceiio et ahl 
l2004ll . 



accretion rate substantially altered the total central lumi- 
nosity, and the decreased central luminosity made the dust 
sublimation radius shrink, that is, a large amount of dust 
condensed in the inner disc region during the fading and 
the subsequent low-brightness state. The bright-state lumi- 
nosity of Lace + Lt ~ 100 Lq implies J?subi ~ 0.7 AU (cf. 
fig. 7 of iDullemond fc Moniiierll2010l ). which might have 
shrunk to ^ 0.4 AU due to a factor of two drop in Lace - The 
observed drop in the ifg-band flux between 2004 and 2007 
suggests similar decrease in the dust sublimation radius, if 
we assume that most of this flux is emitted from the dust 
sublimation zone. Assuming normal interstellar gas-to-dust 
ratio, an excess extinction of Ay ^ 3 mag requires a gas 
column density of ~ 6 x lO'^^ cm~^, which, taking into ac- 
count the size of the dust condensation region, corresponds 
to a volume density of some lO^cm"^. Considering typi- 
cal midplane gas densities an d scaleheights at the dust rim 
(jPullemond fc Monnieill2010l '). the inclination of 62° and the 
fact that the disc of PV Cep is at least 10 times mor e mas- 
sive than a typical HAe star disc (|Hamidouchell201ol ) . such 
densities along the line of sight are likely, and thus the newly 
formed dusty region may account for most of the observed 
optical and near-infrared fading of the star. The observed 
drop in the mid- and far-infrared fluxes indicates that parts 
of the optically thin outer flared disc atmosphere got into 
shadow. 

The transient peak in 2008, according to the colour- 
magnitude diagrams, resulted from variable obscuration 
along the line of sight. A very similar transient peak both 
in amplitude and time-scale appeared in the light curve of 
PV Cep in 1979, f ollowing the end of the outburst (fig. 2 of 
ICohen et al.lll98ll ) , which suggests that the temporary dust- 
clearing might have been related to the end of the outburst. 
A short-lived outflow, similar to that observed in V1647 Ori 
after the end of its outburst in 2006, and probably associ- 
ated with the rearrangement of the ste llar magnetic field in 
response to the dropped accretion rate (|Brittain et al.|[2010l ) 
might have blown out some dust from the line of sight. 

The slowly descending light curve accompanied by red- 
dening between 2008 June and 2009 August may result from 
the dust condensation process. The duration of the low- 
brightness period suggests that the dust, emerging in the 
line of sight, was not confined to dense clumps or warps in 
the disc, but rather it was associated with a restructuring 
of the inner disc due to the changing central luminosity. 

The brightening accompanied by reddening m 2009- 
2010 may indicate that the dust, condensed during the pre- 
vious years, started evaporating. The increased Call fluxes 
suggest that invigorated accretion may play role in the pro- 
cesses. The time-scales of the large-amplitude photometric 
fluctuations suggest that the inner dust rim is involved. No- 
tably, two prominent peaks of the Jc light curve (around 
JD 2,455,103 and JD 2,455,143) are separated by 40 days, 
corresponding to the Kepler-period at a radius of 0.4 AU. 
The low Jfs-band flux, rising fluxes at 3.6 and 4.5 yLxn, as 
well as the reappearance of the nebula suggest that the in- 
ner dust rim of the disc became less efficient in shadowing 
the outer parts of the circumstellar matter. 

This scenario can qualitatively explain the fading in the 
optical, near-, mid-, and far-infrared wavelength regions, and 
the simultaneous variations of the nebula and the Call emis- 
sion line fiuxes. Our results suggest that PV Cep differs from 



both the EXor type stars, characterized by outbursts driven 
by strongly enhanced accretion, and UXor stars, whose vari- 
ability is caused by orbiting circumstellar dust structures. 
The high disc mass and persistent high accretion rate sug- 
gest that PV Cep is a protostar-like object: most of its lu- 
minosity originates from accretion, and its large-scale pho- 
tometric variations reflect the changing inner disc structure, 
resulted from unsteady accretion. Variations in the accre- 
tion rate may be triggered by gravitational instabilities in 
the massive disc of PV Cep. 
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ONLINE SUPPORTING MATERIAL 

APPENDIX A: SPITZER DATA REDUCTION 

IRAC. The IRAC observations of PV Cep were performed 
in two different observing modes. The observations in 2004 
and those in the framework of the post-Helium monitor- 
ing programme were performed using the subarray readout 
mode. The observations in 2006 were performed in high dy- 
namic range readout mode with 0.6 s and 12 s integration 
time. 

The subarray mode allows the utilization of exposure 
time as short as 0.02 s enabling the observation of bright 
sources that would saturate the detector using the other 
readout modes. The data processing was started from the 
Spttzer Science Center (SSC) basic calibrated data (BCD) 



produced by the pipeline version S18.7. The BCD image 
cubes of the 64 frames were combined into two-dimensional 
images using the "irac-subcube-coUapse" IDL routine pro- 
vided by the SSC. The observation performed in 2004 was 
made of nine different dither positions and repeated three 
times resulting in 27 images at each wavelength. In the case 
of monitoring observations four dithering steps were utilized, 
leading to four images at each wavelength. We performed 
aperture photometry on the final images at each wavelength 
using a modified version of the IDLPHOT routines. The aper- 
ture radius was set to 3 pixels, the sky background was com- 
puted in an annulus with an inner radius of 3 pixels and a 
width of 4 pixels. In the course of the sky estimates we 
used an iterative sigma-clipping method, where the clipping 
thresh old was set to 3cr. Following the outline of lHora et al.l 
l|200dl we also applied an array-dependent photometric cor- 
rection and a pixel-phase correction to the measured fiux 
densities. An aperture correction was then performed using 
the values published in the IRAC Data Handbook. The in- 
dividual flux density values measured in the different dither 
positions (at each band) were averaged to get the final pho- 
tometry. The measurement error was estimated as the r.m.s. 
of the individual fiux densities. 

In the case of the high dynamic range mode observa- 
tions, we made the photometry on individual saturation cor- 
rected, CBCD, images. Because of the bright target we used 
only those 12 short exposure time images (with integration 
time of 0.6 s) where PV Cep appeared on the frames. In the 
course of photometry we used the same method as in the 
case of subarray mode observations. 

The final uncertainties of the photometry were com- 
puted by adding quadratically the measurement errors and 
an absolute calibration error of 2% (IRAC Data Handbook). 

IRS. The Spitzer spectroscopic observations of PV Cep 
were obtained in spectral mapping mode with the short low 
(5.2-14.5 ^im, R«60-127), short high (9.9-19.5 ^m, R«600), 
and long high (18.7-37.2 /im, R~600) channels. The spectral 
maps consisted of two nod positions and nine positions per- 
pendicular to the slit. We only considered the two central 
positions, as if it were a normal staring measurement. The 
data were processed at the SSC with the SIS. 7.0 pipeline to 
BCD level. We extracted spectra from the two-dimensional 
dispersed BCD images using the Spitzer IRS Custom Ex- 
traction software (SPICE version 2.1.2). In the case of the 
short low channel, we first subtracted the two nod positions 
from each other in order to remove the sky background, then 
extracted a spectrum from a wavelength-dependent, tapered 
aperture around the star. In the case of the high-resolution 
channels, we extracted spectra from the full slit. Since no ex- 
tended emission could be detected in the IRAC images, we 
calibrated the spectra by assuming that PV Cep is a point 
source. In case the star was not well centered in the slit (i.e. 
the offset from the slit center perpendicular to the slit was 
>0.5"), we corrected the spectra for fiux loss using the mea- 
su red IRS beam profil es, following the procedure described 
in iKospal et al.l (|2008t ). After this, there was still a slight 
mismatch between the different channels. The discrepancy 
could be eliminated by simply adding a linear correction to 
the short high channel. 



MIPS photometry. 24 and 70 nm images of PV Cep were 
obtained as part of the Gould Belt Legacy Survey. We 
started from BCD files reduced with the SSC pipeline ver- 
sion S16.1. At 70 /im, we used GeRT (version S14.0 vl.l 
[060415]) to do column mean subtraction and time filtering 
on the BCD files. We used mopex (version 18.1.5) to create 
mosaics at both 24 and 70 fim, with pixel scales of 2.45" and 
4", respectively. Since PV Cep was saturated at both wave- 
lengths, aperture photometry was not possible. In order to 
determine the fiux of PV Cep, we first obtained aperture 
photometry for 7 or 5 non-saturated stars in the field at 24 
and 70 /im, respectively, (an aperture radius of 7", sky an- 
nulus between 40 and 50", and aperture correction of 1.61 
was used at 24 /im, and an aperture radius of 8", sky annu- 
lus between 39 and 65", and aperture correction of 3.70 was 
used at 70 /im) . Then we normalized the psf profiles of these 
stars to 1 Jy and averaged them. The obtained psf profile 
was then fitted to the non-saturated parts of the profile of 
PV Cep and the scaling factor needed for a good fit directly 
gave the fiux of our target. 

MIPS spectroscopy. The low resolution far-IR (55- 
95 /tm, R« 15-25) spectrum of PV Cep was obtained in the 
MIPS spectral energy distribution (SED) mode. Two ob- 
serving cycles were performed, each consisting of six pairs 
of 3 s long on- and off-source exposures. The on-source and 
the background position was separated by 1'. The data re- 
duction of the MIPS SED observation was started with the 
BCD images (pipeline version S18.12) and the MOPEX soft- 
ware was utilized to perform the necessary processing steps 
(combination of data, background removal, application of 
dispersion solution) and to compile the final image with 
pixel size of 9.8". The spectrum was extracted from the sky- 
subtracted final image using a five pixel wide aperture. As a 
final step an aperture correction was a pplied, the aperture 
correction factors were taken from lLu et al.l ^200^ ). 
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Table A2. Results of the photometry of PV Cep 



JD 


Date 


V{AV) 


i?c{ARc) 


Ic (A/c) 


Telescope 


(2454000+) 


yyyymmdd 


mag 


mag 


mag 




—329 


20(1 "^1 027 


14 5S CO Ifi") 


1 S 28 Cn 06") 


11.78 


(0.06) 


RCC 


—327 


ZrUuo -Luzy 


1 08 To 04^1 


1 "^Q ("0 04^ 


11.90 


(0.11) 


RCC 


—326 


ZUUO -LUOU 




1 '^Q {f\ 04^ 


11.98 


(0.06) 


RCC 


—325 


9nn'^i (^'^^ 

ZUUO ±UO -L 


^d (f) 07~i 


1 "^fi /'O 01 ^ 


11.84 


(0.16) 


RCC 


—308 


9nn'^i 1 1 7 

ZUUO -L J- -L I 


'^Q (f) 04^ 




12.34 


(0.05) 


RCC 


I 






1 ^1 Q8 ('0 09^ 


14.60 


(0.02) 


RCC 


26 


onnfiim? 


17 '^4 ('0 0^^ 


1 c; Q8 ("O 02^ 


14.65 


(0.02) 


RCC 


570 








13.60 


(0.04) 


*-lrfiTTiirlt 


589 


^uuouuuz 






14.11 


(0.02) 




594 


9nn8n'^n7 

ZUUOUiJU 1 






14.25 


(0.04) 


OljIllIllUli 


626 


9nnsnfins 

zuuouuuo 






14.99 


(0.02) 


^ Vi TV1 1 ^1 ^" 

Oi-'lillllLl L 


UZ 1 


zuuouuuy 




1 '^Q (f\ 09^ 


15.03 


(0.02) 


Schmidt 


645 


ZrUUOUUZr ( 




1 '^'^ ("0 09^ 


15.12 


(0.03) 


RCC 


646 


9nnsnfi9s 

ZUUOUUZrO 




1 ^^7 ("0 09^ 


15.14 


(0.04) 


RCC 


647 


9nnsnfi9Q 

zuuouuzry 




1 fil ("0 01 ^ 


15.16 


(0.01) 


RCC 


649 


900^0701 

^UUOU ( U J. 




1 R 70 ("0 01 ^ 


15.21 


(0.01) 


RCC 


650 


^UUOU 1 uz 




1 R R1 f'O 0^1 
xu.ux ^^^J.uo^ 


15.20 


(0.02) 


RCC 


675 


9nnj^n797 

^UUOU 1^1 




1 R 74 ("0 04^ 


15.33 


(0.04) 


OClillllUt 


o / o 


9nn8n798 




1 R 89 /'n 0/1^ 


15.43 


(0.04) 


ociimiUL 


693 


9nnsns"i a 




1 R fi7 /'O 01 ^ 

XU.U / ^U.UX J 


15.31 


(0.01) 




l uo 


9^^s^s9(^ 

ZUUoUoZU 


18 "^7 ^"0 0'^^ 


1 R RQ /'n 0*^^ 


15.09 


(0.04) 


n A 9 9 TYi 

Z.Z-lIl 


707 


9nnsns9s 

zuuouozo 






15.15 


(0.02) 




708 


9nnsns9Q 

zuuouozry 


1 8 04 ^0 (^'^^ 


1 R f^O ("0 0*^^ 


15.14 


(0.03) 


CA 2 2m 


709 


zuuouoou 






15.02 


(0.02) 


^ n Tvi 1 'J" 


710 


9nnKnK'^i 

^UUOUOO -L 




1 R "^f^ ("0 09^ 


15.00 


(0.02) 




738 




1 7 80 /'O 09^1 

J- 1 .OU 1 U.U^ I 


1 R '^9 ("O 01 ^ 


15.00 


(0.01) 


RCC 


757 


9008101 7 


1 7 Q4 ('0 0^"^ 


1 R 41 (C\ 0"^^ 

XU.^X lU.UOJ 


15.07 


(0.03) 


RCC 


758 


90081 01 >^ 


1 7 7^ ('0 09"^ 


1 R 9Q To 01 "1 


14.94 


(0.01) 


RCC 


760 


90081090 


1 7 81 (Cl 02^ 


1 R To 01 ^ 

XIJ.OtJ I U.UX J 


15.00 


(0.01) 


RCC 


782 


90081 111 






14.94 


(0.01) 


OdllllitX I 


808 


2nnsi 207 




1 fi 24 Co OA) 


14.84 


(0.05) 




809 


90081 908 




1 R 1 Q ("0 09^ 
xu.xy yu.yj^j 


14.80 


(0.02) 


^ n vvi 1 *H 'J" 


SAO 


900001 08 






15.18 


(0.02) 


Sclimidt 


841 


900Q01 OQ 


1 8 08 /'O 1 0^ 

±0.\J0 \\J.±\J} 


1 R '=14 f'O 07^ 


15.22 


(0.04) 


SrVimiHt 

OClilliiU. I 


883 


900Q0990 




1 R 84 f'O fl^^ 

XU.O^ I U.UU J 


15.41 


(0.06) 




892 


900Q0294 






15.34 


(0.04) 


SrVimiHt 

kJ Ullll 1 1 U. li 


912 


90000*^91 






15.11 


(0.03) 


OUlillllUt 


916 


9O0Q0'^9'^ 




1 R ^^8 /'O 04^ 


15.24 


(0.04) 


OUlillllUt 


947 


9000049"^ 




1 R 88 ("0 01 ^ 


15.38 


(0.01) 


OUllIIliUlj 


949 


900Q0Z197 




XU. 1 1 ^^u.ux y 


15.31 


(0.01) 


Ot-'llllilU. L 


950 


900Q0zl98 




1 R R9 ("0 09^ 


15.24 


(0.01) 


0*-'lilliii_l. li 


951 


900Q0ZL9Q 




1 R 7Q ("0 09^ 
xu. I y ^^u.uz y 


15.33 


(0.01) 


n YYt 1 ^~ 

kjt-'lilliiU. U 


953 


900Q0'^01 




^f\ 7Q (f\ 09^ 
xu. I y yu.yj.^ J 


15.28 


(0.02) 


RCC 


972 


20090520 




16.97 (0.01) 


15.50 


(0.01) 


RCC 


974 


20090522 




16.97 (0.01) 


15.48 


(0.01) 


Schmidt 


976 


20090524 




16.35 (0.03) 


15.19 


(0.02) 


Schmidt 


978 


20090526 




16.85 (0.02) 


15.43 


(0.02) 


Schmidt 


1001 


20090618 




16.91 (0.04) 


15.44 


(0.05) 


RCC 


1022 


20090709 




16.91 (0.04) 


15.45 


(0.03) 


Schmidt 


1027 


20090714 




16.95(0.04) 


15.49 


(0.04) 


Schmidt 


1058 


20090811 




17.23(0.02) 


15.70 


(0.02) 


Schmidt 


1060 


20090816 




17.24(0.01) 


15.70 


(0.02) 


Schmidt 


1061 


20090817 




17.27(0.01) 


15.71 


(0.02) 


Schmidt 


1062 


20090818 




17.37 (0.01) 


15.86 


(0.02) 


Schmidt 


1071 


20090827 




17.39 (0.01) 


15.78 


(0.01) 


RCC 


1072 


20090828 




17.40 (0.01) 


15.77 


(0.01) 


RCC 



Table A2. Results of the photometry of PV Cep (cont.) 



JD 


Date 


Vc(AV-c) 


i?c{Ai?c) 


/(A7) 


Telescope 


(2454000+) 


yyyymmdd 


mag 


mag 


mag 




1077 


20090902 






15.48 


(0.01) 


RCC 


1093 


20090918 




16.95(0.02) 


15.21 


(0.02) 


RCC 


1094 


20090919 




16.93(0.02) 


15.10 


(0.02) 


RCC 


1095 


20090920 




16.87(0.02) 


15.17 


(0.02) 


RCC 


1096 


20090921 




16.91 (0.02) 


15.17 


(0.02) 


RCC 


1101 


20090926 


18.24 (0.06) 


16.59 (0.01) 


14.93 


(0.01) 


RCC 


1103 


20090928 




16.46 (0.02) 


14.73 


(0.01) 


Schmidt 


1113 


20091008 




17.05 (0.03) 


15.35 


(0.03) 


Schmidt 


1114 


20091009 


18.69 (0.04) 


16.95 (0.01) 


15.32 


(0.01) 


CA 2.2-m 


1115 


20091010 




16.96 (0.04) 


15.48 


(0.06) 


Schmidt 


1116 


20091011 


18.70 (0.04) 


17.00(0.03) 


15.33 


(0.03) 


CA 2.2-m 


1119 


20091014 


18.60 (0.01) 


16.83(0.01) 


15.12 


(0.01) 


CA 2.2-m 


1130 


20091025 




16.85(0.02) 


15.08 


(0.02) 


lAC-80 


1137 


20091101 


18.18(0.04) 


16.61 (0.02) 


14.66 


(0.01) 


lAC-80 


1139 


20091103 




16.54(0.02) 


14.59 


(0.02) 


lAC-80 


1140 


20091104 




16.59 (0.02) 


14.73 


(0.02) 


IAC-80 


1142 


20091106 


18.21 (0.04) 


16.39 (0.02) 


14.46 


(0.02) 


IAC-80 


1143 


20091107 


18.08 (0.03) 


16.40 (0.02) 


14.45 


(0.02) 


IAC-80 


1149 


20091113 




16.51 (0.03) 


14.58 


(0.04) 


Schmidt 


1155 


20091119 




16.95(0.04) 


15.33 


(0.04) 


Schmidt 


1156 


20091120 




16.72(0.03) 


15.12 


(0.03) 


Schmidt 


1159 


20091123 




16.95(0.01) 


15.35 


(0.02) 


Schmidt 


1161 


20091125 




16.98(0.01) 


15.32 


(0.01) 


Schmidt 


1212 


20100115 




16.45(0.05) 


14.53 


(0.06) 


Schmidt 


1213 


20100116 




16.37 (0.06) 


14.46 


(0.07) 


Schmidt 


1253 


20100225 




16 27 (0 13^1 


14.56 


(0.12) 


Schmidt 


1259 


20100303 




16..36(0.13) 


14.38 


(0.12) 


Schmidt 


1288 


20100401 




15.55 (0.07) 


13.76 


(0.05) 


Schmidt 


1295 


20100408 




15.87(0.01) 


13.83 


(0.05) 


RCC 


1331 


20100514 




15.50(0.02) 


13.56 


(0.02) 


RCC 


1372 


20100624 


17.70 (0.06) 


15.81 (0.03) 


14.02 


(0.03) 


Schmidt 


1420 


20100812 


17.45 (0.01) 


15.53(0.01) 


13.64 


(0.01) 


RCC 


1435 


20100826 




16.19(0.04) 


14.31 


(0.05) 


Schmidt 


1461 


20100921 




15.58 (0.03) 


13.53 


(0.03) 


RCC 


1463 


20100923 




15.90 (0.05) 


14.06 


(0.04) 


Schmidt 


1464 


20100924 




15.86 (0.05) 


13.99 


(0.06) 


Schmidt 


1477 


20101007 






14.13 


(0.03) 


Schmidt 


1524 


20101123 




16.15(0.04) 


14.45 


(0.04) 


Schmidt 



Table A3. Spitzer photometry for PV Cep. All fluxes are color-corrected and given in Jy. 



Date 






F4 




Fs.T^im 


Fs.O/am 


F24(im 


Fro^im 


20041029 


4.47 


(0.09) 


6.98 


(0.14) 


8.34 (0.17) 


11.85(0.23) 






20061126 


1.84 


(0.04) 


3.99 


(0.08) 


5.45 (0.11) 


8.16(0.17) 






20070226 














27.12(1.04) 


35.22 (2.13) 


20090916 


3.54 


(0.03) 


5.97 


(0.02) 










20090917 


3.16 


(0.03) 


5.41 


(0.03) 










20090918 


3.19 


(0.07) 


5.31 


(0.03) 










20090919 


3.31 


(0.03) 


5.48 


(0.03) 










20090920 


3.22 


(0.02) 


5.34 


(0.02) 










20090921 


3.11 


(0.04) 


5.19 


(0.04) 










20090922 


3.29 


(0.04) 


5.40 


(0.03) 










20100109 


3.27 


(0.07) 


5.53 


(0.03) 










20100110 


3.14 


(0.04) 


5.31 


(0.02) 










20100111 


3.31 


(0.05) 


5.49 


(0.02) 










20100112 


3.34 


(0.08) 


5.67 


(0.02) 










20100113 


3.26 


(0.04) 


5.46 


(0.03) 










20100114 


3.00 


(0.08) 


5.16 


(0.02) 










20100115 


3.05 


(0.02) 


5.23 


(0.02) 











